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, We perform a general parton level analysis for the search of heavy resonant states 

in the production of tt pairs at the LHC with an integrated luminosity of 30 fb _1 . We 
assume the existence of resonances that only couple to quarks and propose kinematic 
. cuts in order to increase the amount of events produced through quark-annihilation. We 

study the interplay between different variables and their impact on the purity of the 
selected sample. We make focus on the longitudinal (/?) and transverse ipr) momentum 
■ of the tt pair, and the scattering angle (9) in the center of mass reference frame. We 

5_i ! observe that j3 is replaced by 6 as a suitable discriminating variable of quark-annihilation 



processes for invariant masses above 1 TeV. Finally, we illustrate the analysis with a 
gluon resonance of 1.5 TeV and show the improvement in the sensitivity of the signal 
when cuts on 9 are imposed. 
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1 Introduction 



With what seems to be the recent discovery of the Higgs boson [1], the high energy physics 
enters into the solitude of a path that no one knows where or how it may found New Physics 
(NP). Although a light Higgs requires either, fine-tuning, or NP at the TeV scale, the ex- 
perimental results do not show any hint of new phenomena insofar. As a matter of fact, 
natural Supersymmetry, which is one of the most popular theories to avoid fine-tuning and to 
simultaneously explain many other phenomena, is being highly constrained by direct searches 
at the LHC [2]. 

From a general theoretical point of view, one may expect to have NP effects in the detailed 
study of the heavier particles properties. In particular, the top quark is the heaviest known 
particle and its properties have not been explored in depth insofar. In fact, last years results 
from Tevatron [3] seemed to point to NP in the top forward-backward asymmetry. However, 
recent results from the top charge asymmetry at LHC [4,5], have contradicted most of the 
NP proposals except for very few models [6] that have survived. 

Regardless of these recent results in top physics, in this article we consider the inclusive 
pp — > ti process where we study the ti invariant mass (mtt) spectrum as a sensitive observable 
to NP resonant phenomena. This is due to a potential enhancement in the coupling of the 
top quark to the NP. In this work, we address the question of how to increase the sensitivity 
of this observable in the case of NP that couples to quarks and not to gluons. The reason 
for this study is that at the LHC most (> 75 % at y/s = 8 TeV and increasing with energy) 
of the ti events have a gluon in the initial state partons and, therefore, a possible NP as the 
stated above would be per-se diluted just because of the initial state partons of the event. 
For the sake of brevity, from this point forward we call gluon-fusion any event with at least 
one gluon in the initial state. 

The goal of this article is to propose kinematic cuts which enrich the quark-annihilation 
fraction (f qq ) of a selected ti sample in order to increase the sensitivity to the stated NP in 
the m^-spectrum. Moreover, in light of the upcoming experimental results, we focus our work 
to the rati ^ 1 TeV region of energy. This range of energies has two special features which 
determine the results in this work: i) we propose to use as one of the variables to discriminate 
quark-annihilation events, the center of mass scattering angle, which is useful for this purpose 
at high mtt', an d ii) the lack of statistics at high m t i requires a special selection criteria that 
with small cuts makes considerable increases in the sensitivity. This latter is in contrast to 
the low energy regime, where strong cuts that yield a very large increase in the sensitivity 
are sought. 

The search for kinematic cuts that increase f qq in a given ti sample has been studied 
in the last years [7-12] mainly to increase the sensitivity in the top charge-asymmetry at 
the LHC and very few in the m^-spectrum [13, 14]. These articles deal mainly with two 
different features in the production mechanisms that allow to distinguish qq production in 
pp — > ti. The first one is that valence quarks inside the colliding protons tend to have larger 
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fraction of the proton momentum than gluons and anti-quarks. Therefore, top pairs produced 
through quark-annihilation tend to be boosted along the pipe line. This characteristic can 
be measured through the kinematic variable 



which ranges from for not boosted pairs to 1 for maximum boosted pairs. Notice for future 
purposes that the origin of this feature is in the proton parton distribution function (PDF). 
The second feature which allows to isolate top pairs produced through quark-annihilation is 
that initial state gluons tend to produce more initial state radiation (ISR) than quarks. This 
is an effect due to a larger numeric factor in the three gluon vertex in the QCD Lagrangian. 
This characteristic can be measured through the transverse momentum of the top pair, 



The larger is px, the more probable is to have a gluon in the initial state. Notice that, in 
this case, the origin of this feature is pure QCD and is independent not only from the PDF, 
but also from the dynamics in the ti production process. The usage of these two variables to 
increase the sensitivity to NP in the m^-spectrum has been exploited in Refs. [11, 14]. 

There is a third variable which may serve to distinguish quark-annihilation processes. 
This is the center of mass scattering angle (8) of the top quark direction relative to the right 
moving parton. This variable has been previously studied in Refs. [13,15,16]. Since, as we 
show in the article, to take profit of this variable requires large m t i, we use it to study the 
upcoming experimental results where the m t i region above 1 TeV will be better analyzed. 
This variable is related to the dynamics of the pp — > ti process and the spin of the particles. 
The key feature in using the angular distribution is that for large m t i the quark-annihilation 
processes have a smooth angular distribution, whereas gluon-fusion processes accumulate 
most of the events in the forward | cos(6*)| rs 1 region. Or, to be more precise, the t-channel 
gg ti amplitude is the one that peaks the production in the forward region for relativistic 
tops. Therefore, the angular distribution may also be understood as a discriminator of s- 
from t-channel contributions, as in dijet resonance searches [17]. 

The aim of this work is to analyze simultaneously cuts in all three variables (/3, px and 
6) to enhance the sensitivity to resonant NP in the m 4 f-spectrum by increasing f qq -. Notice 
that, although all three variables have different origins, they have some degree of indirect 
correlation through the PDF's. In any case -due to their different origin- none of them can 
be expressed as a function of the other two. The price to pay for increasing f qq is to reduce 
the fraction of the original sample (f s ) that is selected for the analysis. This reduction, which 
always yields an increase in the statistical uncertainty, may not worth the selection. The 
critical point up to where it is convenient to cut the sample depends on how f qq and the total 
uncertainty behave as a function of f s . Henceforth, it is crucial at this point also to take into 
account systematic uncertainties in the analysis in order to know if the selection cut is useful 
or not. However, since a full realistic simulation of the systematic uncertainty would require 
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a detailed simulation of the detector, which is beyond the scope of this work, we use a simpler 
model which is enough to lead us to the desired interplay between f qq - , f s and statistic and 
systematic uncertainties in order to decide whether the selection is suitable or not. 

Under this scenario, for the setup of the problem and its solution, it will be enough 
to consider pp — > ti processes up to parton level including ISR. Hadronization, detector 
simulation and reconstruction would be useful only if one could perform a fully realistic 
detector simulation including its systematic uncertainties. This is left for the experimental 
groups in case they consider to follow the guidelines in this work. The purpose of this article 
is to show that exists an interplay between a selection cut on the three variables (3, px and 
6 and the statistic and systematic uncertainties, that yields an optimal cut which enhances 
the sensitivity of the m^-spectrum to NP that couples to quarks. We show along the article 
that this interplay, and therefore also the optimal cut, is strongly dependent on the energy 
of the process and the accumulated luminosity. For instance, we show that for the 2012 LHC 
data, as the energy increases beyond ~ 1 TeV, strong cuts in /3 should be gradually replaced 
by mild cuts in 8. 

This article is divided as follows. In next section we present the analytic results for the 
angular distributions in the Standard Model (SM) for quark-annihilation and gg-fusion. We 
also study angular distribution of NP models which couple to quarks and not to gluons. In 
section 3, we study the parton level variables /3, pr and 9 and their interplay with f qq and / s , 
as well as the relationships of these two with the uncertainties. In section 4, we present an 
example of a resonant NP in ti production and show how the progressive cuts in the studied 
variables leads to the visibility of the resonant bump. Section 5 contains a discussion on the 
results and previous works, and section 6 the conclusions. 

2 Angular distribution in ti production 

In this section we study analytically the leading order ti production through the processes 
qq — > ti and gg — > ti within the SM and for different possible resonant NP models. At high 
energy, events initiated by quarks have different angular distributions than those initiated 
by gluons. This is due to the gg-fusion amplitude where a top is exchanged in the t-channel 
producing ti pairs mainly in the forward region. Therefore, one may impose kinematic cuts 
to the relevant angular variables in order to disentangle contributions from the two different 
production modes. 

We begin analyzing the angular distribution of the two production mechanisms within 
SM. The expressions for the production squared amplitudes are given in the helicity basis in 
the initial parton center of mass frame. A sum over spin and color for the initial and final 
states as well as proper spin and color averaging factors which have not been included are 
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assumed all along this section. For the quark-annihilation initial state, we have 1 [18] 

\M(qq->tt)\ 2 = 8g 4 (l-fi)sm 2 9, (3) 



LL,RR 

"\M(qq^tt)\ 2 = % 4 (l + cos 2 #). (4) 



LR^RL 

where L(R) stands for left(right) helicity of the t or i, 9 is defined as the angle between 
the direction of motion of the top and the initial parton moving to the right, and ft is the 
speed of the top quark. As it is apparent from these equations, the like-helicity (LL and RR) 
qq ti production is suppressed as ft — > 1. This behavior is a consequence of the helicity 
conservation. 

For the gg initial state we have [18] 

J2\M(gg^tt)\ 2 = M + I S (1 ~ ffX* + $ + $ ^ °)> (5) 



LL^RR 

E 

LR,RL 



\M(gg^tt)f = fg* J + ^^ 2 /3^in 2 6>(l + cos 2 g). (6) 



For ft close to 1, the common factor to both equations governs the growth of ti production 
through gg-fusion as |cos#| — > 1. Furthermore, this growth becomes stronger as /3j — > 1. 
Hence, we may expect that the angle # becomes a more useful variable to distinguish quark- 
annihilation from gg events for increasing values of the invariant mass of the ti pair. Again, 
like-helicity contribution is suppressed in the limit (3 t — > 1 as for quark-annihilation pro- 
duction, except for |cos#| > 0.9 where both production mechanisms contribute. This is 
understood because the t-channel that arises in this limit allows to flip the helicity of the 
tops. Notice that at threshold (f3 t = 0) both production mechanisms have an isotropic angular 
distribution. Then, no angular distribution based discrimination between the two production 
mechanisms is expected close to the threshold. 

In Fig. 1 we show the normalized like- (left) and unlike- (center) helicity angular distri- 
butions from Eqs. (3)-(6) for three values of ft = 0.57, 0.94, 0.98 {m tl = 425 GeV, 1050 GeV 
and 1700 GeV, respectively). Also, we plot the total angular distributions for the sum of 
all final helicity states (right). Notice that for quark-annihilation production the normalized 
like- and unlike-helicity angular distribution are independent of (3 t since in Eq.(3) (3 t depen- 
dence enters just as a global factor and in Eq.(4) there is no dependence on (3 t at all. We 
can see that qq — > ti production in an unlike-helicity state is more likely to be produced in 
the forward region whereas qq — > ti production in a like-helicity state is mostly central. The 
suppression in the like-helicity qq — > ti production for increasing values of m t i results in a 
slight population of the forward region as we can see from the total qq — > ti angular distri- 
bution. On the other hand, for j3 t ^ the unlike-helicity production for gg — > ti is turned 



1 As it is known, the study of the squared amplitudes is sufficient for the analysis developed in this section 
since the phase space in the two-body angular differential cross sections is independent of cos 6 itself. 
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Figure 1: [color online] From left to right: normalized unlike-, like-helicity and total angular 
distributions for qq — > ti (top) and gg —> ti (bottom) for three values of f3 t = 0.57 (solid 
blue), 0.94 (dashed red), 0.98 (dot-dashed green). These correspond to m t i = 425 GeV, 1050 
GeV and 1700 GeV, respectively. For quark-annihilation production the normalized like- and 
unlike-helicity angular distribution are independent of (3 t (notice that this only takes place 
when the angular distributions are normalized). 

on and both like- and unlike-helicity production begin to populate the forward region as /3 t 
increases. Hence, although both mechanisms tend to populate those regions for increasing 
values of m t t, gg — > ti production makes it further significantly due to the particular behavior 
of its total angular distribution as fit — > 1. The comparison between qq — > ti and gg — > ti 
productions in this limit is shown in the right-hand side column of Fig. 1. These differences 
between the angular distribution of the two production mechanisms for large m t i will prove 
to be useful to distinguish them. We may expect from these plots that as /3j — > 1 a loose 
upper cut in | cos#| results in a slight decrease of the total number of events together with a 
sudden increment of the quark-annihilation fraction of the sample. This is in opposition to 
f3 where one usually imposes tight cuts that end up removing a large amount of events from 
the sample in order to get an enhancement in the quark-annihilation fraction. 

In the presence of NP that couples only to quarks, the angular distribution of the gg-fusion 
events remains the same as in the SM. Therefore, the fraction of these events that are removed 
with a given cut on 9 does not change with respect to the SM. For instance, we obtain from this 
simple parton level analysis that making a forward angular cut |cos#| < 0.85 for = 0.98, 
only 35 % of the gg- fusion initiated events pass the cut. This has to be compared with 
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Figure 2: [color online] Angular distribution of the qq —> ti NP events through spin-0 (a), 
spin-1 vectorial or axial (b) or spin-2 (c) massive resonances. The different curves correspond 
to different f3 t as stated in Fig. 1. 



the fraction of SM quark- annihilation events that are left after the same selection, which is 
~ 80 %. Moreover, it is easy to verify that this selection cut is also useful for NP propagating 
in the s-channel. We test this with three benchmark cases where the NP corresponds to 
spin-0, spin-1 or spin-2 massive resonances as is shown in Fig. 2. Since the spin-0 colorless 
resonance (Fig. 2a) does not interfere with the quark-annihilation QCD production and has 
no angular dependence, the same cut yields that 85 % of the NP quark-annihilation events 
are selected. For the spin-1 resonance (Fig. 2b) we choose a color octet G' with vectorial 
couplings. The amplitude squared in this case is 

J2\M(qq g/G' -+ ti)\ 2 = A(s)(2 - ft + ft cos 2 9). (7) 

spins 

The same angular dependence that the SM is obtained with A(s) being a factor that depends 
on the resonance mass and width. This is consequence of the equal angular dependence of 
the QCD and NP amplitudes. Thus, we also obtain that 80 % of this NP quark- annihilation 
events are left after the selection. Moreover, we find no difference with this and the case 
of an axial spin-1 color octet. (Notice that in this case, if we had chosen 9 as the angle 
between the top and the quark instead of a fixed beam direction, as is usual in forward- 
backward asymmetry studies, we would have obtained an asymmetric angular distribution 
for the qq — > G' — > ti process, but the symmetric cut |cos6*| < 0.85 would have had the same 
effect.) Finally, in Fig. 2c we illustrate as an example a spin-2 colorless resonance [19] which 
we restrict its couplings to quarks. In this case we find that 74 % of the NP events pass the 
selection. Thus, provided that the NP propagates in the s-channel, the cut in 9 could be 
used to improve the quark- annihilation fraction and the sensitivity to resonant NP. In fact, 
although we find some slight differences in the efficiency of a given cut on |cos#| to keep 
quark-annihilation events from NP, we always obtain the same amelioration in reducing the 
dilution coming from gg-fusion events. 
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3 Using /3, Pt and 6 to enhance the sensitivity to NP 
resonances 



We now discuss results from simulations of ti events within the SM in order to find phase 
space regions where quark-annihilation is favored over gluon-fusion as the dominant produc- 
tion mechanism. We consider each kinematic variable separately to analyze them as potential 
filters of quark-annihilation processes. Finally, we perform random scans of cuts on these vari- 
ables and show that to impose restrictions on 9 is more efficient to isolate quark-annihilation 
production than to force limits on /3 for the high m t i region. 

We have simulated SM inclusive ti production at the LHC@8TeV up to one extra jet with 
MadGraph5 [20] and then showered it with Pythia [21]. To avoid double counting, we have 
matched the matrix element to the parton shower through the MLM scheme [22] implemented 
in MadGraph5. This generation procedure allows to include nontrivial distributions. It is 
beyond the scope of this work to deal with final states, therefore we have assumed a simplified 
selection cut \r]t,t\ < 2.5 (for future purposes, it is worth noticing that rj = 2.5 corresponds to 
cos(#) = 0.987, for the case where both laboratory and center of mass frames coincide). The 
prediction for the inclusive production cross section has been obtained with MCFM LHC@8TeV 
at NLO as 225.2 pb [23] and an overall selection efficiency for semileptonic ti detection of 
6 % [5] has been used in the analysis. This strategy has ended up with a total of 400.000 ti 
events at the 201 2's 30 fb" 1 . 

We show in Fig. 3 the differential distribution of quark-annihilation (red dashed) and 
gluon-fusion (black continuum) for each one of the kinematic variables /3, px and 9 for three 
different bins of m t t (400-450 GeV, 1000-1100 GeV and 1600-1800 GeV) that will be used from 
now on all along the analysis. Notice also that the central values of these bins correspond to 
the m t t values selected in the previous section as referential cases. Since all the distributions 
have been normalized to one, the plots are not useful to compare the number of events 
initiated by quark-annihilation or gluon-fusion. Instead, they allow to evaluate whether a 
given cut is suitable for increasing the quark-annihilation fraction. Thus, we can see from the 
first row of the figure that quark-annihilation events tend to populate the high j3 region more 
than gluon-fusion events. In particular, the peak in the distribution of quark-annihilation 
events at large /3 in the 400-450 GeV can be understood according to the following. For 
large mtt, due to the suppression of the PDF's, it is easier the energy to be shared by the 
initial quark and anti-quark than having a configuration with large momentum difference 
between them. This situation is relaxed as m t i decreases and it is more likely to find a 
quark with a large momentum. We conclude from the plots in the first row that a lower 
cut in f3 always increases f q q. By contrast, from the second row, we expect that an upper 
cut in px may favor the quark-annihilation production. Finally, the angular distributions 
for both production mechanisms is depicted in the third row. As we have discussed in the 
previous section, both quark-annihilation and gluon-fusion events have essentially the same 
almost isotropic distribution in the 400-450 GeV bin. For increasing m t i values, the angular 
distributions of both mechanisms starts to differentiate from each other. For the 1000-1100 
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M = 400 — 450 GeV M = 1000 - 1100 GeV M = 1600 - 1800 GeV 

(94200 events) (2200 events) (120 events) 




Figure 3: [color online] Differential tt distribution for quark-annihilation (red dashed) and 
gluon-fusion (black solid) events for each one of the kinematic variables f3, pt and 6 in three 
invariant mass bins (400-450 GeV, 1000-1100 GeV and 1600-1800 GeV). All the distributions 
have been normalized to one. 

GeV bin and even more remarkably for the 1600-1800 GeV bin, gluon-fusion events produce 
t and t around the incoming direction of the initial partons whereas the events initiated by 
quark-annihilation do not have such a strong 9 dependence. 

We notice that the efficiency in performing cuts on and px to disentangle quark- 
annihilation from gluon-fusion events is slightly sensitive to m t i. On the contrary, from 
the angular distributions we observe that gluon-fusion events populate the forward direction 
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for increasing m t i values and 9 becomes a more convenient variable to discriminate quark- 
annihilation events. In order to go into greater detail on this matter, we have plotted in Fig. 4 
the quark-annihilation fraction (f qq ) vs. the fraction of events (/ s ) that remain after a given 
cut for each variable. For the 1600-1800 GeV bin, a larger quark-annihilation fraction with 
a less loss of events in the sample is achieved by setting cuts on 9. Therefore, it is expected 
that cuts on 9 to be more suitable to small samples of events as it is the case in resonance 
searches beyond 1 TeV. In summary, cuts on 9 may lead to larger values of f s in the relevant 
invariant mass region. 



M = 400 - 450 GeV M = 1000 - 1100 GeV M = 1600 - 1800 GeV 




Figure 4: [color online] f qq (vertical axes) vs. f s as the sample is cut in (3 > (3 C (red dotted), 
Pt < PTc (blue dashed) or |cos#| < cos# c (magenta solid) for three invariant mass bins 
(400-450 GeV, 1000-1100 GeV and 1600-1800 GeV). The range of the cuts are f3 c G [0,0.95], 
pre E [20,300] GeV and cos9 c G [0.1, 1] 

So far we have considered each kinematic variable separately. Now, we combine cuts in 
the three variables with the aim to improve the quark-annihilation fraction. We perform 
then simultaneously a lower random cut (3 > (3 C and upper random cuts px < Pre and 
| cos#| < cos9 c . We have set f3 c G [0,0.9], p Tc G [20,300] GeV and cos9 c G [0.1, 1]. A cut on 
f3 stronger than 0.9 may suffer from large systematic uncertainties [10] and p Tc = 20 GeV [5] 
corresponds to the minimal experimental sensitivity. 

In Fig. 5 we show f qq vs. f s for a scan on the three variables (vertical lines area), and only 
in (3 and px (shaded area) (namely, no cut on 9 is demanded). We have plotted for the three 
invariant mass bins 400-450 GeV, 1000-1100 GeV and 1600-1800 GeV. As expected from the 
discussion given in the previous section, an enhancement in f qq is obtained at large invariant 
masses when we combine cuts on all the three variables instead of only imposing cuts on /3 
and pt- Moreover, we have also checked that, for the 1600-1800 GeV bin, the maximum f qq 
for a fixed f s is reached without any cut on (3. The independence on the cut on f3 is consistent 
with the behavior observed in Fig. 4 within the same m t t range. Therefore, in order to end 
up with f s as large as possible for a given value of f qq , we conclude that performing cuts 
on 9 is more efficient than restricting (3. This is an important observation since f3 is usually 
thought as a central kinematic variable in phenomenological analyses that aim at achieving 
quark-annihilation and gluon-fusion discrimination [7-12]. 

These SM results can be used to search for new physics in the it invariant mass spectrum. 
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For the 1600-1800 GeV bin, simulated NP distributions follow the same behavior as in right 
plots in Fig. 3. In particular, to impose upper cuts on | cos#| is still an efficient procedure 
to remove a considerable fraction of events initiated by gluon-fusion while keeping a large 
f s , as discussed in previous section. Therefore, the results of this section provide us with a 
guideline to look for new physics in the m t f-spectrum. 



M = 400 - 450 GeV M = 1000 - 1100 GeV M = 1600 - 1800 GeV 




fs fs fs 



Figure 5: [color online] f qq (vertical axes) vs. f s for a scan on the cuts on j3, pr and 9 (vertical 
lines area), and on j3 and px (shaded area) (no cut on 9) for three invariant mass bins (400-450 
GeV, 1000-1100 GeV and 1600-1800 GeV). 

We consider now a possible matter of concern since, as the value of f qq becomes larger, 
the number of events in the sample decreases with a logical increment in the statistical 
uncertainty. This may spoil the potential power of the whole analysis and, for resonance 
searches, we do not only need a large f q q but also to have uncertainties under control. From 
Fig. 5 we can see that, as f s gets reduced, there is always a set of cuts that improves f qq . 
If the systematic uncertainties dominate over the statistical ones, any cut enhancing f qq 
always improves the sensitivity to new physics. However, if we are searching for resonances 
in the region beyond 1 TeV, a reduction on the amount of events is expected and further 
cuts may lead immediately to a regime where both systematic and statistical uncertainties 
are competitive. Even if f qq is large, the statistical uncertainties can spoil the sensitivity to 
NP [14]. Therefore, a regime where systematic and statistical uncertainties are of the same 
order has to be considered as a minimal necessary condition to be fulfilled. 

It is beyond the scope of this work to perform an exhaustive analysis of the uncertainties 
and, for the sake of the discussion, we simply assume that the relative systematic uncertainty 
is constant. This is a good approximation as far as we do not push the cuts to extreme values. 
We model the systematic uncertainty in each bin as a syst = cN, where N is the number of 
events in the bin and c is the relative systematic error in corresponding bin. A typical 
value for c in a tt -spectrum measurement without jet substructure top tagging is about 
20 % (c = 0.2) [24]. (This number, and the selection efficiency, are expected to change if top 
tagging through jet substructure [25] is used in the analysis.) Thus, we get that the systematic 
and statistical uncertainties are of the same order provided that N ~ 1/c 2 ~ O(10). It is 
worth noticing at this point that if the systematic uncertainties are reduced (c decreases) 
then the optimal number of events for a given bin increases in order to reduce the statistical 
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uncertainties. Therefore, if the systematic uncertainties are reduced then the optimal selection 
cuts should be relaxed. 

We illustrate the above discussion with an example. At 8 TeV with 30 fb _1 , we expect in 
the SM 120 events in the bin of 1600-1800 GeV of invariant mass. As we have mentioned, the 
systematic and statistical uncertainties will be comparable for N ~ 25, i.e. f s ~ 0.2. From 
Fig. 5 we observe that for f s « 0.2 there exists a set of cuts on (3 (which turns out to be 
needless, (3 > 0), px and 9 that achieves f qq ~ 0.75. Whereas for cuts only on (3 and px we 
would need f s « 0.1 to achieve the same f qq (see the dashed line in Fig. 5c). Accordingly, we 
would obtain the same significance with a cut on (3, Pt and 9 at 30 fb _1 as with a cut only on 
f3 and pt for 60 fb _1 . For this reference example, the set of cuts is (3 > 0, pt < 24 GeV and 
|cos#| < 0.89 for the former case, and (3 > 0.38 and px < 20 GeV for the latter. This case 
shows that to include 9 as a discrimination variable is equivalent to double the luminosity. 



4 Application of the studied cuts in a specific example 

In this section we apply the previously discussed cuts to a specific NP model in order to 
explicitly show the amelioration in the invariant mass spectrum sensitivity to a new resonance. 
We study the effects of including cuts in different sets of /3, pt and in 9 to visualize the 
improvement in sensitivity as more variables are included in the cuts. 

We take as an example a benchmark model with a gluon resonance of 1.5 TeV reminiscent 
of strongly coupled/composite models. A characteristic feature of this kind of models is that 
this massive gluon couples stronger to the heavy than to the light quarks, and it cannot be 
created from a fusion of two QCD gluons. Thus, this kind of models can illustrate and test 
the studied cuts. The interaction with SM particles is: 

n 

where q stands for all quarks. We take f( u ,d, s ,c) L R = 0.025, f tR = 4.6 and all other couplings 
are set to zero. With this choice of parameters the width of G' yields Tq' = 290 GeV, and 
the resonance would not be visible with the LHC running at 8 TeV and an accumulated 
luminosity of 30 fb~ x without any cut to improve the f qq fraction. 

With the same computational tools used in the previous section, we have simulated a 
30 fb" 1 sample of pp — > ti and assumed a reference selection efficiency of 6 % [5]. We have 
simulated only SM and jointly SM+NP, to take into account the interference coming from the 
large width effects. In Fig. 6 we show the comparison in the region of interest of the expected 
mff-spectrum of only SM (histogram with error bars, in black) and SM+NP (histogram, in 
red) for: (a) no cuts at all; (b) cuts in (3 and px\ (c) cuts in /3, px and 9; and (d) an ideal cut 
where only the quark-annihilation events in the sample are left. The error bars correspond to 
statistical and systematic uncertainties, as discussed in the previous section. We choose the 
cuts in plots (b) and (c) from a scanning that looks for the cuts that maximize the significance 
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Figure 6: [color online] it invariant mass spectrum for the 30 fb _1 expected LHC sample. 
The (black) histogram with error bars corresponds to only SM, and the (red) histogram 
corresponds to SM+NP. In panel (a) we plot both spectrum without applying any cut on the 
sample, in (b) we apply the best cut that can be achieved in using (3 and p?, in (c) the same 
as in (b) but using (3, px and 9, and finally (d) corresponds to an ideal cut, as explained in 
the text. In each panel we determine the significance of the deviation according to a x 2 test 
(see text). 



of a x 2 test to 3 bins of 100 GeV width from 1500 GeV to 1800 GeV. We have checked that 
if we had chosen the a priori cuts from the previous section SM analysis to the 1600 — 1800 
GeV bin, then we would have obtained an analogous, but not optimal, improvement. 

To quantify the gradual improvement in resolving the resonance in Fig. 6, we show in 
Table 1 the f s , f q q , p- values for the x 2 tests, and the Gaussian-equivalent significance (a) of 
all the set of cuts in the figure. We also show the value of the cuts for Fig. 6b and Fig. 6c. 
As we can determine from the table (second and third row), a cut in /3 is not too tight at 
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plot 


cuts 


/. 


fqq 


p- value 


a 


(a) 


w/o cuts 


1.00 


0.41 


0.36 


0.4 


(6) 


/3 > 0.45, p T < 75 GeV 


0.20 


0.63 


0.13 


1.1 


(c) 


(3 > 0.08, p T < 65, |cos#| < 0.66 GeV 


0.22 


0.75 


0.02 


2.1 


(d) 


ideal cut 


0.41 


1.00 


0.004 


2.7 



Table 1: [color online] Set of cuts used to differentiate SM+NP from SM using different set of 
variables. The first two columns indicate the plot in Fig. 6 and the performed cuts, the third 
and fourth columns show their effect on the sample (f s and f qq ), and the following columns 
the expected differentiation in computing a y 2 test, as explained in the text. 



these energies and luminosities. Moreover, in passing from second to third row, when the cut 
in 9 is included, we see that the cut in (3 is practically replaced by a cut in 9. We see in this 
example that this exchange of cuts increases f s by ~ 10 % and f qq by ~ 20 %, which yields 
an improvement in the p- value from 0.13 to 0.02. 

We have repeated the same computations as in this section for other models of resonant 
NP and obtained similar results. For the best cuts in the analogous to Fig. 6c we always 
obtain (3 > — 0.1, | cos(6 l )| < 0.65 — 0.85 and p T < 30 — 70 GeV. This is in agreement with 
the previous section analysis. It could be valuable to stress at this point that the results 
in this section have the only purpose of showing the qualitative effect of these cuts in the 
increment in the sensitivity through a specific example. 



5 Discussions 



The aim of this section is to relate the present analysis with previous works. Besides px and 
j3, the rapidity of the top and spin correlations of its decays also have been used to separate 
NP and SM qq — > ti production from the gg — > ti background. First, we discuss in this section 
some differences among the present analysis that makes focus on 9 as a suitable kinematic 
variable and previous studies where the rapidity y is used to impose tight selection cuts. And 
finally, we discuss briefly about the different regimes where 9 cuts and spin correlation are 
useful to improve the f qq fraction. 

It is well known that there is a correlation between the center of mass angle 9 and the 
difference of rapidity Ay = y l — y l , 

1 + &cos0 
1 — p t cos 9 



where f3 t = \/l — 4m^/s. The partonic energy s entering this relation may dilute the corre- 
lation between this two variables. Only if /3 t is fixed, there is a one-to-one correspondence 
between 9 and Ay. However, if we look at a bin of invariant mass where /3t has large varia- 
tions then there could be a lower correlation between them. For instance, if we consider the 



14 



m t i bin of 1600 to 1800 GeV, Ay varies up to 4 % when 9 is fixed. Therefore, for a high m t i 
bin, cuts on 9 and Ay are roughly equivalent. Moreover, if we take several bins in the high 
m t i region, for the same cut on 9, we will have a different (but equivalent) cut on Ay in each 
bin. On the other hand, close to the threshold, flt varies greatly and the correlation is more 
diluted. Thus, within the 400-425 GeV bin, Ay varies up to 20 % for fixed values of 9. 

Cuts on | Ay | have already been considered in the literature [9,10,26]. However, they 
have been applied at low energies and for certain class of models. For instance, in Z' flavor 
violating models, it are produced through the exchange of the Z' in the t-channel and then it 
produces it events mostly in the forward region. This situation, which does not correspond 
to resonant NP, is quite the opposite to the present analysis. For such models, it has been 
shown that a lower cut on \Ay\ enhances the sensitivity to the new interactions. Since these 
interactions arise at low m#, there is not a considerable increment of the gluon-fusion fraction 
that may spoil the enhancement generated by the new interactions. In cases where the new 
physics contributions are significant at higher m t t, we have to deal with an interplay between 
increasing the NP t-channel contribution and decreasing the sensitivity with the increment 
of the gluon-fusion it production. 

The rapidity has also been used in the context of the charge asymmetry in it production 
at the LHC [8,11]. In those cases, the basic idea is also to discriminate quark-annihilation 
from gg-fusion events by imposing different cuts on y l and y l . The procedure traces the 
discrimination back to the correlation between y and /3. However, the angular variable 9 (or 
Ay) plays no role in that method. 

Finally, for many of the models that have been proposed to accommodate the Aps, there 
exist a high m f j resonance that induces an interference at low m t t, for instance, a heavy 
axigluon that has been integrated out [27]. Therefore, the new physics contribution arise 
at low mtt an d this may seem to be the reason why 9 has never been taken as a relevant 
kinematic variable in previous studies. 

Spin correlations in the production and decay of tops have been extensively discussed [9,18, 
28]. The aim of those studies is to compare predictions made within the SM or its extensions 
with the measured angular distributions of the decay products of t and i. In general, different 
models give place to distinct angular distribution patterns and the measurement of spin 
correlations may provide additional and relevant information on the structure of possible 
new interactions. Such spin correlations are optimized when the asymmetry between like- 
and unlike-helicities reaches certain maximum value. Usually, for a given axis from which the 
spin projection is defined, one seeks to perform kinematic cuts that increase this asymmetry. 
This in turn allows to discriminate like- and unlike-helicities. 

From Eqs. (3)- (6) and their corresponding discussion in Section 2 we have learned that 
angular distribution for both quark-annihilation and gg-fusion production mechanisms are 
essentially flat close to the threshold. Therefore, 9 cuts cannot be used to improve the quark- 
annihilation fraction. On the contrary, a distinction based on spin correlation is envisaged 
since only like-helicity gg — > it production is possible whereas the two spin production mech- 
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anism are present for qq — > ti as we can see from the Eqs. (3)-(6). The opposite situation we 
encounter at large m t i . There, only unlike-helicity qq — > ti and gg — > ti production mech- 
anisms are possible, being hardly expected that spin correlation being useful to distinguish 
quark-annihilation from gg-fusion production contrary to 9 cuts that has been proven to be 
suitable to increase the quark-annihilation fraction. 

6 Conclusions 

We have studied the collider phenomenology of ti resonance searches above 1 TeV at 30 fb _1 . 
We have focused our analysis on resonances that couple exclusively to quarks which, at the 
LHC, would be naturally hidden behind the large gluon-fusion ti production. 

The usual variables to enlarge the quark-annihilation fraction f qq - of a ti sample are the 
longitudinal and transverse momentum of the ti pair, (3 and pt, respectively. In addition 
to these variables, we exploit the fact that for masses above 1 TeV the large boost of the 
tops in the center of mass system yields a correlation between their center of mass scattered 
angle 9 and the initial partons of the event (this correlation is suppressed for invariant masses 
below 1 TeV). Since the origin of the relationship between the initial state partons and (3 is 
the proton PDF, whereas for is the initial state radiation, and for 9 the dynamics of the 
partonic collision, there are not a priori reasons to expect these variables to be considerably 
correlated, and all three should be used together to texture the initial state partons of an 
event. 

We have simulated ti production in the large invariant mass regime and scanned over 
different set of cuts on the variables mentioned above to increase the quark-annihilation 
fraction of the sample. Given a set of invariant mass bins where the NP resonance could be 
expected, the optimal cut is the one that enhances the most the significance of a resonance 
search over the SM background. This is achieved when a set of cuts increases f qq while keeping 
the selected fraction f s still large since, for resonances beyond 1 TeV, a critical reduction in 
the amount of events could spoil the signal due to an increase of the statistical uncertainties. 

Our results show that j3 is not as useful as expected in this resonance search. The reason 
for this is that this variable requires tight cuts in order to yield significant increases in f qq , 
and there are not too many events in the large invariant mass bins in the 2012 data. On 
the other hand, 9 does not reach large values of f qq , but it has an important increase as one 
begins cutting the sample. Hence, the reasons why 9 is useful in the large invariant mass bins 
are two-folded: the variable requires relativistic tops to distinguish the initial state partons, 
and in this region the cuts should be moderate. At last, the cut in p T is always useful for 
any invariant mass bin. 

We have studied a reference case of a 1.5 TeV resonance and found that the optimal 
cuts in these three parton level variables is in agreement with the previous sections analyses. 
We have also shown how the significance of the resonance search behaves as one gradually 
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includes all variables in the selection cuts. Summarizing, when we include all three variables, 
we found that /3 is not relevant, and px and 9 provide the most useful cuts to enhance the 
significance of the resonance search. 

The results in this article are a guideline to point the searches of NP resonances above 1 
TeV for 30 fb _1 at a center of mass energy of 8 TeV. It is important to stress the perspectives 
for upcoming luminosity and energy increases. For NP resonances coupled to quarks and 
masses above 1 TeV, a rise in the luminosity will make 9 remain as a useful variable, whereas 
f3 will improve as a quark-annihilation discriminator too. On the other hand, a rise of the 
center of mass energy, will allow the exploration of heavier resonances within this framework. 
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